Lattice theory of trapping reactions with mobile species 
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We present a stochastic lattice theory describing the kinetic behavior of trapping reactions 
A + B ^ B, in which both the A and B particles perform an independent stochastic motion 
on a regular hypercubic lattice. Upon an encounter of an A particle with any of the B 
particles, A is annihilated with a finite probability; finite reaction rate is taken into account 
by introducing a set of two-state random variables - "gates", imposed on each B particle, 
such that an open (closed) gate corresponds to a reactive (passive) state. We evaluate here a 
formal expression describing the time evolution of the A particle survival probability, which 
generalizes our previous results. We prove that for quite a general class of random motion of 
the species involved in the reaction process, for infinite or finite number of traps, and for any 
time t, the A particle survival probability is always larger in case when A stays immobile, 
than in situations when it moves. 

PACS numbers: 05.40.-a, 02.50Ey, 82.20.-w 



I. INTRODUCTION 
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Kinetics of chemical reactions involving diffusive 
species have attracted a great deal of scientific interest 
since the pioneering work by Smoluchowski [1]. Since 
then, many novel and conceptually important results 
have been obtained [2-6]. In particular, it has been 
proved in specific cases that the classical, mean field 
chemical kinetics does not apply, at least in low dimen- 
sional systems [4-11]. 

Trapping A + B B reactions (TR), involving ran- 
domly moving A and B particles which react " when they 
meet" at a certain distance R, provide an example of 
chemical reactions showing a pronounced deviation from 
the text-book predictions. 

For the TR two situations were most thoroughly stud- 
ied: the case when As diffuse while Bs are static, and the 
situation in which the As are immobile while Bs diffuse 
- the so-called target annihilation problem (TAP). In the 
case of static, randomly placed (with mean density b) 
traps the A particle survival probability Pa(^) shows a 
non-trivial, fluctuation-induced behavior [8,9,12-17] 
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which is intimately related to many fundamenal problems 
of statistical physics [8,9,12-18]. 

Survival probability Ptarget{t) of an immobile target 
A of radius R in presence of point-like diffusive traps B 
(TAP) can be calculated exactly for any d (see Refs. [21] 
and [4,8,19,20]): 
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where D = Db and Ks{t) is the d-dimensional 
Smoluchowski- type " constant" , defined as the flux of dif- 
fusive particles through the surface of an immobile sphere 
of radius R. Decay forms in systems with hard-core in- 
teractions between Bs [22] or with fluctuating chemical 
activity [23] have also been discussed. 

On contrary, the physically most important case of TR 
when both As and Bs diffuse was not solved exactly. It 
has been proven [9] that here PA[t) obeys 
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which equation defines its time-dependence exactly. On 
the other hand, the factor Xci{Da, Db) remained as yet 
an unknown function of the particles' diffusivities and 
d. Since the time-dependence of the function on the rhs 
of Eq.(4) follows precisely the behavior of drKgir), 
one might expect that the SA provides quite an accurate 
description for this situation and following its spirit to set 
Da — supposing that traps diffuse with the diffusion 
coefficient D = Db + Da- As a matter of fact, it has been 
often tacitly assumed that when both of species diffuse 
PA{t) obeys Eq.(2) with (t)^^\t) defined by Eq.(3) and 
D = Da+Db- On the other hand, it has been shown that 
Xd{DA, Db) is less than the corresponding prefactor in 
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Ks{t) [31] and that it may be bounded by a non- analytic 
function of Da and Db [24]. A perturbative approach for 
calculation of Xd{DA,DB), as well as corrections to the 
SA in ID systems were presented [20]. It has been also 
noticed that Xci{Da, Db) is not a function oi D = Da + 
Db only, since the diffusion-reaction equation are not 
separable [20] . This lack of knowledge of the precise form 
of Xd{DA, Db), of course, constitutes an annoying gap in 
the general understanding of the fluctuation phenomena 
in chemical kinetics. 

Recently, some very interesting and unexpected results 
have been established for trapping A + B ^ B reactions 
involving randomly moving species [25,26], which have re- 
solved, at least in part, this problem. It has been shown 
that in one or two dimensions [25], or more generally in 
systems, in which the fractal dimension of the B par- 
ticle trajectories is greater than the dimension d of the 
embedding space [26], (i.e. in case of the so-called "com- 
pact exploration" [27]), the leading at long times kinetic 
behavior of perfect trapping is essentially independent of 
the A particle diffusion coefficient. In other words, it has 
been shown [25] that in such low dimensions, the leading 
long-time decay of the A particle survival probability in 
systems in which the A particle diffuses and the decay in 
systems in which it is fixed at the origin are exactly the 
same. 

The derivation of this rather surprising result relies 
heavily on the assumption that the A particle has a larger 
probability to survive until a given time t (at least when 
t — * oo) if it stays immobile rather than when it moves 
randomly. In Ref. [25] some arguments have been pro- 
posed in favor of this conjecture, based on the analysis 
of the decay exponents in systems with a finite number 
of Bs. Subsequently, it was shown rigorously in Ref. [28] 
that it is indeed the case in a one-dimensional continuum; 
on the other hand, the claim that Ref. [28] presents a rig- 
orous proof of this conjecture for d = 2 docs not seem to 
be justified; as a matter of fact, Eq.(2) of Ref. [28] does 
not make sense for d = 2 and the logarithmic correction 
does not follow from it, since particle's radius is not taken 
explicitly into account within the approach of Ref. [28] . 

On the other hand, in a recent article [29], this con- 
jecture has been examined rigorously within the con- 
text of reactions between the particles executing ran- 
dom walks on c?-dimensional lattices; here, such a con- 
jecture has been referred to as the "Pascal principle", 
since it is reminiscent of a famous philosophical asser- 
tion of Blaise Pascal, who claimed that "all misfortune 
of man comes from the fact that he does not stay peace- 
fully in his room" [30]. In Ref. [29], we showed that as 
t —>■ oo the Pascal principle-like inequality between the 
survival probabilities of a diffusive and of an immobile 
A particles is valid in any dimension, provided that the 
A particle performs some rather general continuous-time 
jump process on a hypercubic lattice, while the B parti- 
cles perform independently a discrete time lattice jump 



process, which also satisfies some rather natural assiimp- 
tions. The same conclusion was also obtained in Ref. 
[29] for a much more general case of stochastically-gated 
reactions, which mimic situations with finite elementary 
reaction act constants. 

We also emphasize that very similar Pascal principle- 
like inequality has been proven earlier for the process of 
hopping transport of an excitation on a disordered ar- 
ray of immobile donor centers in presence of randomly 
placed, immobile quenchers [31]. We note, as well, that 
recent results obtained for the ballistic A+A ^ annihi- 
lation process [32] are compatible with such a principle. 

On the other hand, the analysis in Ref. [29] is rather 
condensed and moreover, some of the assumptions in- 
voked, as well as some of the constraints imposed on 
particles' random walks, seem to be unnecessary and thus 
can be safely relaxed. Consequently, our purpose here is 
to complete the proof of the Pascal principle-like inequal- 
ity between the survival probabilities of the diffusive and 
immobile A particles and to extend it in several direc- 
tions. In particular, we proceed to show that the Pascal 
principle-like inequality holds at any finite time t, as well 
as for both infinite and finite number of traps. More- 
over, we shall consider here the case when the chemical 
activity of the B particles fluctuates in time between ac- 
tive and inactive states. We set out to show that the 
Pascal principle also applies for this much more complex 
and realistic situation. As in our previous work [29], we 
will focus here solely on the lattice formulation of the 
model. The continuous-space case, which requires much 
more delicate analysis, will be studied elsewhere. 

The paper is structured as follows: In Section 2 we for- 
mulate the model, introduce basic notations and define 
the properties of reaction and random walks executed by 
the species involved. In section 3 we will focus on the 
reaction kinetics in case of perfect trapping; that is, on 
the case when the A particle gets annihilated with proba- 
bility 1 upon the first encounter with any of the traps B. 
Further on, in section 4, the Pascal principle-like inequal- 
ity between the survival probabilities of a diffusive and 
immobile A particles will be extended to the case when 
the chemical activity of the B particle fiuctuates ran- 
domly between active and inactive states, which mimics 
more realistic situations in which an annihilation of an 
A particle upon its encounter with any of the Bs takes 
place with a finite probability. We consider here a rather 
general case when such activity fluctuations can be cor- 
related in time. Some intermediate calculations, as well 
as analysis of the behavior in some special cases are rele- 
gated to the Appendices A and B, in which, in particular, 
the special case of Polya random walks is considered. 
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II. LATTICE MODEL OF TRAPPING 
REACTIONS BETWEEN MOBILE SPECIES 

It is well-known that lattice models of diffusion- 
controlled reactions yield, at least for sufficiently large 
times, the kinetic laws that are essentially the same as 
those obtained within the continuous-space descriptions. 
Thus in the present work we shall consider a lattice model 
of trapping reactions, which will simplify significantly 
our analysis. One of advantages of such a consideration, 
apart of the fact that it allows for much more lucid anal- 
ysis than in the continuous-space limit, is that we are not 
forced to attribute to particles a finite, non-zero radius, 
which allows to consider the behavior in systems of any 
spatial dimension. Finally, for sake of simplicity, we will 
restrict our analysis here to hypcrcubic lattices; most of 
the results, of course, could be readily extended to other 
types of embedding lattices. 

Consider Na particles of type A and Nb particles of 
type B, which are initially placed at random at the sites 
of a rf-dimensional hypercubic lattice, containing M sites. 
All particles perform independent jump processes on the 
nodes of the lattice. Each A particle can be destroyed (in 
the general case with a finite probability) as soon as this 
A particle appears on the same lattice site simultaneously 
with any of the B particles. The B particle remains 
unchanged after the reaction event, which corresponds 
to the annihilation mechanism 

A + B^B, (5) 

and represents the customary trapping-like reaction. One 
can also envisage a more general catalytic reaction pro- 
cess of the form: 

A + B^C + B (5') 

where the catalyst B promotes the transformation of an 
A into some product molecule C, the product molecule 
C being immediately extracted from the system. 

In regard to the reaction probability, we will distin- 
guish between two situations: the one of perfect trapping, 
or purely diffusion-controlled trapping, in which case any 
A gets annihilated with probability 1 upon the first en- 
counter with any of Bs, and that of imperfect trapping 
for which the annihilation of the A by any B takes place 
with a finite probability < 1. To mimic this condition, 
we will introduce a set of additional random variables, 
attached to each B particle, which will describe their 
instantaneous reactive activity. Finally, we will assume 
in what follows that collisions (simultaneous encounters) 
between two (or more) A particles are possible and do 
not affect these A particles or their random walks, and 
similarly, that collisions between the B particles are pos- 
sible and do not lead to any reactions. In other words, 
neither A nor B particles have hard-core interactions and 
no single-species reactions may take place. 



Now, as far as particle motions are concerned, we face 
here the following problem: on one hand, in regard to 
dynamics of A and B particles, we have to define two 
different random processes with different characteristics, 
e.g. diffusion coefficients, which may be used afterwards 
as tunable parameters. On the other hand, these ran- 
dom processes must allow a rigorous analysis, which is 
not always the case. If we choose, for example, that 
both species perform random hopping motion in discrete 
time, then it will be quite difficult to work out a rigorous 
formalism in which two random processes have different 
"waiting" times at lattice sites. If, on contrary, we choose 
that both processes evolve in continuous time, then we 
will face purely mathematical difficulties in treatment of 
the events in which the particles of different species ap- 
pear simultaneously at the same lattice site. Not for- 
saking the generality, we thus choose here the "mixed" 
case, which seems to us most suitable for the rigorous 
description. 

We thus assume that all A particles perform identi- 
cal and independent continuous-time jump processes, so 

that they can jump at any time moment from one lattice 
site to any other site. No other hypothesis or additional 
constraints on their motions are required. 

Further on, we suppose that all Bs perform identical 
and independent discrete-time random walks; that is, at 
any integer time n > 0, any B particle can jump from a 
lattice site y to site y' with a given probability p{y'\y), 
where y' can be identical to y, i.e. B can remain at the 
site which it occupies at time moment n. Let denote 
the position of a given B particle at time moment n > 0. 
We assume that this random walk satisfies the following 
conditions: 

(i) the random walk is homogeneous in space and in 
time: the probability p(2/'|y) is independent of time and 
of the initial position, such that 

p{y'\y) =p{y' -y)- (6) 

(ii) the uniform distribution p{y) = 1/M is stationary 
for the B particles. We assume that it holds at time 

0, and hence, it is realized at all times. This condition 
implies the bilateral normalization relation 

y 

We remark that a stronger condition would be to assume 
that the probability p(j/'|j/) satisfies the detailed balance, 

1. e. 

p{y'\y) =p{y\y') 

which, together with the condition in eq.(2), implies that 

p{y-y')=p{y' -y)- (8) 

However, the condition in eq.(3), which clearly follows 
from eq.(4), is sufficient for obtaining our main results. 
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(Hi) at any time n, the conditional probability P{Yn = 
y\YQ = 0) of finding a given B particle at an arbitrary 
position y at time moment n, provided that it started its 
random walk at the origin, docs not exceed the return 
probability P{Y„ = OlYo = 0). That is, 

P{Yn = y\Yo = 0) < P{Yn = 0\Yo = 0)=Rn (9) 

This condition seems to be quite plausible for any sym- 
metric random motion in a uniform medium if eq.(4) 
holds, since here P{Yn = y\yo = 0) is invariant upon 
reversal y into —y. Hence, this probability should always 
have an extremum for y = 0, which is likely to be a max- 
imum. It should be noticed, however, that inequality in 
cq.(5) docs not hold exactly for the usual Polya random 
walk, when the particle jumps at one of neighboring sites 
at each integer time moment. For instance, on a one di- 
mensional lattice of unit spacing, the inequality in cq. (5) 
is not satisfied if n and y are odd, since in this case the re- 
turn probability i?„ = 0. Nevertheless, one easily obtains 
(without any significant lack of generality) random walks 
satisiying condition (Hi): for instance, one may consider 
a modified Polya random walk on a rf-dimensionnal hy- 
percubic lattice, such that at each integer time moment 
a walker has a probability po to remain at the site it 
occupies, and a probability (1 — po)/2d to jump at one 
of the neighboring sites. In this case, it can be shown 
that the inequality in eq.(3) is verified if po > 1/2. More 
generally, eq.(5) holds for any homogeneous and sym- 
metric probability j)(a; I y) if p(x|.t) = p(0) > 1/2 (see Ap- 
pendix A). Furthermore, our conclusions can be extended 
to cover the case of the Polya random walks, as shown in 
Appendix B, where this special question is discussed in 
detail. 

We finally remark that in the continuous-space case, 
in which the random walk is replaced by a Brownian 
motion, such a question does not arise at all, since here 
the probability density is always maximal and centered 
aroimd the initial position. Thus the inequality in cq.(5) 
appears very naturally in unbiased diffusion problems, 
but it can also be verified for non symmetric jump prob- 
abilities. 

We close this section by adopting some conventions on 
how to introduce reaction events into the model. We 
assume that a given B particle can only annihilate A at 
integer times n > 0. If at a non integer time A jumps 
on a site which is occupied by a particle B, it will be 
only annihilated at the next integer time n. Finally, we 
remark that the probability that an A particle performs 
a jump exactly at an integer time is 0, which allows to 
neglect consideration of such events. Note also that all 
these assumptions won't change the global behavior of 
the system. They thus merely serve for convenience of 
exposition. 



III. PERFECT TRAPPING ON A LATTICE. 

We consider first the case of perfect trapping in which 
case an annihilation of an A particle; takes place at the 
first encounter with any of B particles. Our aim here is to 
demonstrate, in a rigorous way, the Pascal principle-like 
assertion that the survival probability of an A particle 
which moves randomly on a lattice is less or equal to the 
survival probability of an immobile A particle. 

A. Mean-field kinetics of the trapping reaction. 

We start with a reminder on the predictions of a con- 
ventional mean-field approach [2]. One notices first that, 
clearly, the average number < N^in) > oi A particles 
surviving up to an integer time n is the sum of proba- 
bilities that a given A particle survives up to this time 
moment n. Since all of them have identical evolution 
laws, one has 

{NA{n)) = iVA(0)*(n), 

where 'l'(n) denotes the survival probability of a single 
particle A. Since the particles B are completely insensi- 
tive (as far as their motions are concerned) to particles 
A, ^{n) can be evaluated independently for each particle 
A. Thus, it is legitimate to consider only the survival of 
a single A particle in presence of N particles B. 

In terms of the conventional mean-field kinetics [2] , one 
obtains then an exponential decay form for \E'(n): 

*(n) = exp {-khn), (10) 

which should hold in any dimension d. In the last equa- 
tion k is the reaction constant and b stands for the mean 
density of the B particles. 

Note that in case of perfect trapping eq.(lO) becomes 
senseless, since here k = oo. Indeed, it has been well- 
known for a long time, both for the continuous-space and 
lattice models, that the decay law in cq.(lO) does not 
hold, at least for d = 1 and d — 2 [1,4,7-11], so that 
mean-field approach fails and a detailed stochastic theory 
is needed. 

To illustrate the deviations from the mean-field behav- 
ior in eq.(lO) and the actual decay forms, let us consider 
the case when one of the species only is moving [7,8]. In 
this illustration, we follow closely the methods outlined 
in Ref. [33]. 

B. Survival probability of an A pEirticle. 

Let us call Ta the A particle trajectory, and first 
suppose that it is given. Then, we denote as .tq = 
0,xi, . . . ,Xn the successive positions of the A particle at 
the integer times to = 0,ti, . . . ,tn- 
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We suppose next that the waiting time of A at each 
lattice site, i.e. the time which an A particle spends on 
this site between successive hops, is a stochastic variable, 
so that two successive positions arc not necessarily dif- 
ferent, and not necessarily nearest neighbors. 

Further on, we denote the ith B particle, i = 1, . . . ,N, 
as Bi and as - the stochastic trajectory of this par- 
ticle. Next, let Qi{n\FA) be the conditional probability 
that Bi does not destroy A up to time moment n, for a 
given trajectory Ta of A. Because all BiS move and act 
independently of each other, the conditional probability 
^'(n|rA) that the particle A survives up to time moment 
n for a given Ta, factorizes 

N 

^{n\TA) = l[Qi{n\rA) (11) 

i=l 

and hence, the overall A particle survival probability 
obeys 

*(n) = (*(n|rA))r„ (12) 

the average being taken over all possible trajectories of 
A from t — to t = n. Furthermore, since all B particles 
are identical, one has that Qi(n|r^) = Q{n\TA) for all i, 
and hence 

*(n) = {Q{n\rAf)r. (13) 

where, once again, the average is being taken over all 
possible trajectories F^ of the A particle. 

C. The survival probability in the thermodynamic 
limit. 

Let us denote Yo,Yi, . . . ,Yn the successive positions of 
a given B particle at time moments 0, 1, . . . , n, and Fj^^ 
- a trajectory starting from Iq = 2/0 at time 0. One can 
write then 

Q{n\TA) = {Q{n\TA,yo))yo, (14) 

where Q(n\TA,yo) stands for the conditional probability 
that a given B particle, starting its random walk from 
position yo at time moment 0, does not destroy A until 
time moment n for a given trajectory F^. The brackets 
< . . . >y(, in eq.(14) denote averaging with respect to all 
possible initial positions yo of a given B particle. 

We now assume that the probability of the initial po- 
sition Yo is uniformly distributed among the M available 
sites. Then we have 

{Q{n\rA,yo))yo = -^^QHTA,yo) 
yo 

and eq.(9) can be written 



*(r^) = ({1 - ^ E (1 - yo))}^)r. (15) 

Vo 

Turning next to the thermodynamic limit, i.e. setting 
N —>■ oo and M ^ cxd, while keeping their ratio fixed, 
N/M — > 6, 6 being the concentration of the B particles, 
one obtains for the A particle survival probability at time 
n the following expression 

*(«) = (exp{-6E(l-Q(n|FA,yo))})r. (16) 

yo 

Hence, the survival probability \l/(n) is simply related to 
the probability that a given B, starting from yo, destroys 
A at some time t < n, for a given trajectory of A, which 

is 

P{n\TA,yo) = l-Q{n\TA,yo) (17) 

Similar results were obtained [4,7-11,26,33] in the par- 
ticular case when the A particle is immobile, i.e. for the 
so-called target annihilation problem. In this particu- 
lar case there is no averaging over F^ as in the previous 
formulas, and the integral reaction rate is thus defined 

by 

K{n\rA)=J2P{n\TA,yo), (18) 
yo 

which replaces in this case the term kn of the conven- 
tional kinetic law in eq. (6) . On contrary, in more realistic 
situations when A also moves, the average over the tra- 
jectories Ta makes the explicit calculation of the survival 
probability impossible in most cases. 

D. A basic inequality. 

Let us define P^{k\TA,yo) as the conditional probabil- 
ity that B, starting from j/o at time 0, meets A for the 
first tim,e at time A;, given the trajectory F^. Then, the 
conditional probability P(n|F^, yo) that B, starting from 
yo, destroys A at or before time moment n is given by 

P(n|r^,yo)= P\k\TA,y<,) (19) 

0<fe<n 

The conditional probability that the trajectory of B (ex- 
tended after the possible annihilation of A) meets F^ at 
time n (not necessarily for the first time) satisfies the 
equation 

P(y„ = Xn\Yo = yo) = 

P{Yn=Xn\Yk = Xk)P\k\TA,yo), (20) 

0<fe<n 

where P{Yn = Xn\Yk = Xk) = ^x^.x*, and P^(p\TA,yo) = 
0. 
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Summing both sides of the last equation over all ini- 
tial positions yo and using the relation in eq.(3), which 
applies to P{Yn = a;„|lo = Vo), we obtain 

1= P{Yn=Xn\Yk = Xk)S{k\rA), (21) 

0<k<n 

where we have used the notation 

S{n\TA)=Y.P\n\TA,yo) (22) 

Vo 

Next, using the inequality in eq.(5) we obtain from 
eqs.(15) and (18), the following basic inequality 

1< Rn-kSik\TA), (23) 

0<k<n 

where Rn-k is the probability of return to the starting 
point in n — fc steps, which is a well-known quantity for 
all classical random walks. 

We notice that if the A particle stays immobile, the in- 
equality in eq.(19) becomes the equality, since here Xk = 
for all times k. 

Now, let F{s) denote the generating function of some 
function F{n), 

F(5)=^i^(n)s" (24) 

n>0 

Multiplying both sides of the inequahty in eq.(19) by s" 
and performing summations, we have then 

(y^) < Ris)S{s\rA), (25) 

where R{s) is the generating function of the return prob- 
ability Rn, while 

§(s|r^) = ^ (^pi(n|r^,2/o))s" (26) 

n>0 \ 3/0 / 

Note that again, the inequality in eq.(21) becomes the 
equality in the particular case when A is immobile, so 
that 

^(s|0) < SislTA) (27) 

where S{n\0) denotes the S{n\rA) in case when A is im- 
mobile. 

On the other hand, one readily notices from eqs.(12) 
to (14) and (18), that the A survival probability at time 
n is just 

*(n) = (exp(-6K(n|r^)))r^ (28) 

where 

K{n\TA)= Y S{k\TA) (29) 

0<fe<n 



In the limit n ^ oo this expression coincides formally 
with the generating function of S, if s — > 1, which sug- 
gests that inequality (23) corresponds, at least asymp- 
totically, to the similar inequality 

K{n\0) < K{n\TA) (30) 

In this inequality the right-hand side corresponds to the 
case of an immobile particle A. Consequently, the in- 
equality in eq.(26) implies that the annihilation is faster 
if A moves than if it is immobile, in agreement with the 
Pascal principle. 

However, the generating functions "^(slO) and S'(s|rA) 
tend to 00 when s — > 1, and the derivation of eq.(26) 
requires a more careful analysis, which is the purpose of 
the next paragraph. 

E. General form of Pascal principle. 

Let us turn back to the inequality in eq.(19) and recall 
that it becomes an equality in the case when the A par- 
ticle does not move. Then, we may formally rewrite the 
inequality in eq.(19) in the following form: 

0< Y Rn-k[Sik\TA)-S{km (31) 

0<fe<n 

Next, let us introduce two auxiliary functions L„ and 
M„, such that 

L„ = S{k\TA) - S{k\0) (32) 

and 

M„= ^ Lk = K{n\rA)-K{n\0) (33) 

0<fe<n 

By definition, we have Rq = 1 and Lq = 0. Then, the 
inequality in eq.(31) can be straightforwardly written as 

Mn> Yl {Rn-l-k - Rn-k) Mk (34) 
l<fe<n-l 

Now, it can be shown (sec Appendix A) that Rk is a de- 
creasing function of k. Asuming that it has been proved 
that Mfc > for all < A; < n - 1, it follows from eq.(34) 
that Mn > 0, and inequality in eq.(31) is proved by in- 
duction, since Mq = 0. 

Consequently, for any time n, the A particle survival 
probability ^'(n), defined by eq.(24), in case when the 
A particle does not move is less or equal to the survival 
probability in case when A is mobile. This result is much 
stronger than the asymptotic form of the Pascal principle 
obtained in Ref. [29]. 

We note also that this result holds for any finite sys- 
tem with a finite number of traps (i.e. not necessarily in 
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the thermodynamic hmit). Here, the survival probabil- 
ity is given by eq.(9), in which equation Q{n\TA) denotes 
the probability that a given B particle does not meet A 
before or at time moment n. In fact, the probability 
P^(n|r^) that a given B particle meets A for the first 
time at time moment n obeys 

P\n\TA)=Q{n-l\rA)-Qin\TA) (35) 

and we have 

Q(n|r^) = l- J2 P\k\TA) (36) 

0<fe<n 

In Section §3.3. we introduced P^(fc|r^) as the con- 
ditional probability that B, starting from yo at time 0, 
meets A for the first time at time k. Particles B arc uni- 
formly distributed among the M lattice sites at time 
(and at all time as well), so that, with the notations used 
in eq.(18), we find that 

P\k\TA) = ^^pi(A;|r^,2/o) = ^S{n\rA), (37) 
yo 

while eq.(35) reads 

Q{n\TA) = 1 - ^^("ir^) (38) 

Thus, the inequality in eq.(26) implies that, whatever 

may be the number of B particles, particle A has a 
higher probability to escape if it stays immobile, than 
if it moves. Note that this conclusion had been drawn 
previously by Bray and Blythc [25] for systems with a 
finite number of traps within the context of survival of a 
mobile pray A in presence of a finite number of predators 
B. Within this context, the Pascal principle- like inequal- 
ity in eq.(26) appears to be even more sound. 

To close this section we note that the inequality in 
eq.(26) may be questioned for usual Polya random walks 
[34], since the condition in eq.(5) is not strictly fulfilled. 
It is shown in Appendix B how our results can be ex- 
tended to this case. 

IV. IMPERFECT TRAPPING. 

A. Time correlated chemical reactivity fluctuations 

We now modify the model presented in Section 2, as- 
suming that the A particle has a finite probability (which 
may depends on time) to survive when encountered by 
a B particle. This case occurs if the reaction is not 
purely controlled by diffusion: at each encounter, an- 
other stochastic process arises and allows the reaction to 
be eventually completed, or to fail. This process is an 
elementary reaction act. 



If there is a single A particle, it is physically plausi- 
ble to assume that at each of its encounters with any 
of the Bs, the latter can be either in a passive internal 
state with a (possibly time-dependent) probability p{t) 
((0 < p{t) < 1)), or in an active state with probability 
1 — p{t). In the latter case, the A is destroyed, whereas 
it remains intact if B is passive and they may harmlessly 
coexist until the B changes its reactive state. In Ref. [29] 
we have already addressed this problem, assuming that 
this reaction probability was constant and independent 
of all prior events. In many circumstances, however, this 
assumption is not justified, and, in particular, the sur- 
vival probability of the A particle during its encounter 
with any of the Bs may itself depend on the trajectory 
of these particles. We will not treat this difficult problem 
in general, but only consider a special situation, in which 
the A particle survival probability depends on some inter- 
nal, stochastic property of particle B. Then it is possible 
to take into account the time correlations of its fluctua- 
tions. This model can be justified as an approximation of 
certain phenomena, such as possible fluctuations in the 
activity of the catalyst in a chemical reaction [36]. 

More precisely, we assume [33,37] that each particle 
B can be in a passive state or in an activated state 
1, the waiting time Tj in state i being a stochastic vari- 
able independent of prior events, distributed following an 
exponential law 

P{Ti >t)= exp (-A,t) (i=0 or 1), (39) 

where Aq and Ai are given positive constants. 

Now, the transition probability for the internal state 
I(t) of B is then given by the well-known "random tele- 
graph" law [35] 

P{I{t) = i\m = j) = ai + (% - ai)e-^* (40) 

with 

A = Ao + Ai and ao = Ai/A, ai = Aq/A (41) 

Thus, ao = p is the asymptotic probability that A 

survives a collision with B, whereas ai = 1 — p is the 
asymptotic reaction probability at each encounter. The 
elementary reaction act constant k, mentioned in the be- 
ginning of Section 2, is just k ^ (1 — p)/p. We logically 
assume that the internal state of B is initially in its sta- 
tionary probability distribution, as well as at the first 
encounter of with A, but at the next encounter the law 
given by eq.(40) should be used. 

Extending eq.(16) to the present situation with a 
stochastic elementary reaction act, we find 

arP{Y,,=Xn\Yo = yo) = P}{n\TA,yo) + 
+ («i+aoe-^("-'=^) X 

0<k<n-l 

X P{Yn = Xn\Yk = Xk)P\k\rA, Vo) (42) 
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In fact, aiP{Yn ~ a;„|yo = yo) is the probability that 
B meets A at time n while it is in its active state: the 
probability for B to be in its active state at time n is inde- 
pendent of the trajectories of B or A, and is given by the 
stationary value o.\ , since no value of the internal state 
is assigned before time n. Furthermore, P}{k\TA,yo) is 
the probability that B meets A in its active state at time 
k for the first time after 0, with P}{0\TA,yo) =0. If -B 
meets A in its active state at time n, then necessarily the 
same situation occurred for the first time at some time 
k, < k < n. If k < n, then the probability for B to be 
again in its active state is given by eq.(40), which gives 
rise to the last term in the right hand side of eq.(42). 

Equation (42) can be written in a more compact form 



aiP{Yr, = Xn\Yo = yo)= ^ (ai + aoe-^("-'=)) 

0<k<n 

X P(y„ = Xn\Yk = Xk)P\k\rA, Vo) 



(43) 



Summing both sides of it over the initial position y^, we 
then obtain 



ai = ^ (ai +aoe-^("-'=)) 

0<fe<ri 

P{Yn = Xr,\Yu = Xk)S{k\^A) 



(44) 



where S{k\T a) is the probability that the annihilation of 
A occurs at time fc, for a given trajectory of A. 

Using next the inequality in eq.(5), we find the follow- 
ing relation 



1< y (l + %-A(n-fe) 



Rn-kS{k\TA) (45) 



We now make use of the generating functions technique 
and obtain, returning to the notation p = ao, that 



1 



1 



< 



R{s) 



P 



1-p 



R{se-^) 



S{s\Ta) 



(46) 



which becomes the equality in case A is immobile. The 
conclusions follow as previously: the generating function 
of the conditional reaction probability at time n is mini- 
mal if A is immobile, i.e., 



s{s\o) < s{s\rA) 



(47) 



Consequently, the integral reaction rate K_{n\TA) is min- 
imal if A stays immobile 



K{n\0) < K{n\TA), 



(48) 



the left hand sides of eqs.(47) and (48) denoting the quan- 
tities corresponding to an immobile A, respectively. 

The fact that eq.(48) holds at any time n can be proved 
directly by using the inequality in eq.(46) exactly in the 
same fashion as it has been done in Section 3.3. (when 
all collisions are reactive (p = 0)). 



B. Asymptotic reaction kinetics 

Let us first consider the special case when A is immo- 
bile. Then, the asymptotic kinetic behavior follows from 
eq.(46): 



(1 - s)S{s\0) = 



P 



R{s) + -^R{se-^) 
1 — p 



(49) 



In one and two dimensions, R{s) tends to infinity when 
s ^ 1, so that the terms due to the reactivity fluctuations 
in the right hand side of eq.(49) do not affect the kinetics, 
which proceeds exactly in the same way as in the case of 
perfect trapping reactions. 

In three dimensions R{s) tends to a finite limit 1/S 
when s — > 1, where 5* is the probability that a given B 
particle never returns to its initial position (see Appendix 
A). Then, the left-hand-side of eq.(49) tends to an effec- 
tive, "apparent" reaction constant k, which satisfies the 
inverse addition relation 



P 



1 _ 1 

k~S^ 1-p 



Rie-') 



(50) 



if A is immobile [33,37]. 

Equation (50) shows that if A is immobile, the reac- 
tion rate k is an increasing function of the relaxation fre- 
quency A of the activity fluctuations, so that the survival 
probability decreases with A, if p is maintained constant. 
It can be shown that this remarkable property is more 
general and is also valid if both particles are mobile [37]. 
In the case of an infinite relaxation frequency, or uncor- 
related fluctuations, eq.(50) becomes 



1 _ 1 p 



(51) 



Equations (50) and (51) are particular cases of the 
"inverse addition law" which is well-known in chemistry 
[11,36,39]. In fact, such a law is valid if the reaction can 
be considered as a succession of independent steps, which 
is the case for uncorrected fluctuations. It was discussed 
in this context in our previous works [36]. 

We shall now partially extend these results for a mobile 
particle A. In fact, in one or two dimensions it has been 
shown [25,26] that when A is annihilated as soon as it 
meets any of B particles (p = 0), the survival probability 
does not depend on the motion of A in the limit n — > oo 

*(n;p = 0) w ^o{n;p = 0) 

where is the survival probability in case of an immo- 
bile A. 

If A has a finite probability to survive at each en- 
counter, the overall survival probability ^'(n;p. A) is ob- 
viously larger than if p = 0. However, we have shown in 
eq.(47) that '^{n]p, A) is smaller than the survival proba- 
bility \E'o(n;p, A) in case of an immobile particle, and the 
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latter is asymptotically independent of p. Thus, we can 
write 

'^o{n;p= 0) « = 0) < *(n;p, A) 

< *o(n;p,A) «*o(n;p = 0) (52) 

so that for large n 

* (n; p, A) « *o (n; p, A) « *o (n; P = 0) (53) 

which shows that in one and two dimensions the reac- 
tivity fluctuations of B as well as the motion of A do 
not affect the survival probability of A in the asymptotic 
limit n ^ oo, except if B is immobile, in which case 
the survival probability has a very different and unusual 
behavior [13]. 

In three dimensions, on contrary, the fluctuations actu- 
ally change the reaction kinetics. The survival probabil- 
ity decreases exponentially and the overall reaction rate k 
is given by the inverse addition law in eq.(50). However, 
it is unclear if such a law still holds when both parti- 
cles move, since then we have only proved inequality in 
eq.(53). Thus, the first and the last approximate equal- 
ities in eq.(51) do not hold in three dimensions, whereas 
the relations in eq.(51) are valid. It is known [9] that 
for large n, ^{n;p = 0) decreases exponentially, as well 
as the survival probability of an immobile particle A, 
^'o(n;p, A), so that it may be assumed that ^(n;p. A) 
also decreases exponentially with a constant larger than 
k given by eq.(50), but it is difficult to estimate this con- 
stant precisely. 

V. CONCLUSION 

We have developed the stochastic lattice theory of the 
annihilation kinetics of a species A by another species B, 
in systems in which the A and B particles perform in- 
dependent, stochastic motions which can be rather gen- 
eral. We obtained formal expressions for the survival 
probability of A. This probability cannot be evaluated 
exactly if A is actually mobile. However, we proved that 
Thu A particle survival probability is always larger in 
case when A is immobile than when it moves. We have 
shown that this the so-called Pascal principle holds for a 
large class of stochastic motions, provided B executes a 
random walk satisfying certain reasonable assumptions. 
This conclusion is of a special importance in view of its 
implications on chemical reactions or population dynam- 
ics. It also allows to demonstrate that in low dimensions 
the survival probability of the A particle is essentially in- 
sensitive to its motion and fiuctuations of the reactivity, 
and does not obey the conventional mean-field laws of 
chemical kinetics. This result is extended to the case of 
stochastically-gated reactions, including the case when 
the fluctuations of reactivity are time-correlated. Fur- 
thermore, the method used here allows to obtain the 



chemical constant of a stochastically-gated annihilation 
of immobile A particles in a straightforward manner. 

The stochastic analysis of chemical reactions should 
be developed in different directions, in order to consider 
more realistic models. In particular, it would be interest- 
ing to address the case when the activity of A can also 
fluctuate. However, the most necessary improvement of 
the theory would be to extend it to the analysis of anal- 
ogous reaction kinetics in continuous space and time. 

APPENDIX A: PROPERTIES OF THE RANDOM 
WALK OF PARTICLES B 

Stochastic motion of B particles. The constraints 
imposed on this stochastic process are described in Sec- 
tion 2.1. We show here that the main condition in eq.(5) 

P(y„ = y\Ya = 0) < P(y„ = OlFo = 0) = i?„ (Al) 

is satisfied, if the elementary transition probability is 
symmetric and obeys 

p{x\y) = p{x — y) = p(y — x) and p(0) > 1/2. 

To show this, it is expedient to use first the well-known 
formula for the propagator of a random walk on a d- 
dimensional regular lattice (see, e.g. Ref. [34]): 

P(y„ = y\Yo = {)) = -!-[ e-^-y{Hz)rdz, (A2) 
{'^^) Jb 

where B is the first Brillouin zone of the lattice, while 
<f){z) is the so-called structure function 

<^(^) = ^e->(t/) (A3) 
y 

in which equation z ■ y stands for the scalar product of 
two d-dimensional vectors z and y, while dz represents 
the differential element in a d-dimensional space. 

Now, according to our assumption, p{y) = p{—y), 
which implies that 

?^(^) = 51 ■ y)p(y) = p(^) + 

y 

+(1 - p(0)) ^ cos {z ■ y)p{y)/{l - p(0)) (A4) 

Evidently, the second term in eq.(A.4) is bounded from 
above by unity if p(0) > 1 — p(0), or p(0) > 1/2. In this 
case, P(F„ = y\Yo = 0) = (27r)'' cos (^ • y)(.^(z))"dz 
is maximal for y = 0, which proves the inequality in 
eq.(A.l). 

Furthermore, one readily notices that if this condition 
is fulfilled, Rn is a decreasing function of n. 

Relation with first return time. The probability 
R]^ that the first return of B to its initial position occurs 
at time moment n is classically obtained from the relation 
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P{Yn = 0\Yo = 0) = R„ = Rl+ Rn-kRl = 

l<fe<n-l 

= Rn-kRl forn > 1 (A5) 

0<k<n-l 

where Rq = 1 and Rq = 0. Then, the generating func- 
tions of Rn and R!^ satisfy 

R{s) - 1 = R{s)R^{s), (A6) 

which yields 

R{s)^{l~R\s))-^ (A7) 

S being the probabihty that the B particle never returns 
to its initial position. 

APPENDIX B: THE CASE OF POLYA RANDOM 
WALKS 

The assumptions of section 2 exclude the Polya random 
walks, or any random walk such that there is a prob- 
ability to stay immobile at each integer time: p(0) = 0, 
i.e. a random walk in which a particle is forced to make 
a move at each integer time moment. In this case, the 
probability to return to the initial position is obviously 
at any odd time moment, and the inequality in cq.(5), 
which plays a basic role in our reasonings, holds only at 
even moments of time. More precisely, possible displace- 
ments of the random walker in this case can be divided 
into two complementary subsets Eq and Ei, such that 
the total displacement during time n necessarily belongs 
to £^0 if n is even, and to Ei if n is odd. Thus Pascal 
principle cannot apply in a strict sense. 

We can recover the Pascal principle for Polya random 
walks if we slightly change the rules of our model, impos- 
ing, for instance, that the A particle moves only on the 
lattice Eq, and that the B particles are distributed on Eq 
only. Thus no reaction can occur at odd times, and we 
only consider even times n = 2n' . Then, the inequality 
in eq.(5) applies, as well as all previous calculations, and 
the Pascal principle holds. 

However, it is interesting to discuss the case when B 
performs a Polya random walk, if the evolution of parti- 
cles is not restricted on a sub-lattice. The eq.(16), 

P{Yn = Xn\Yo = yo) = P{Yn = Xn\Yk = Xk) 

0<k<n 

P\k\TA,yo) (Bl) 

is still valid, (with possibly many vanishing terms), but 
the inequality in eq.(19) cannot be deduced from it. 

If, in eq.(B.l), n — fc is even, Xn — Xk should belong to 
Eq, and inequality in eq.(5) holds 

P{Yn = Xn\Yk = Xk) < Rn-k (B2) 



On the contrary, if n — is odd, Rn-k = 0, but we 
have 

PiYn = Xn\Yk = Xk) = Ypi^ri\y)P{Yn-l = ^l^fc = Xk) 

y 

V 

Consequently, from eq.(B.l) we can deduce the in- 
equality 

P{Yr, = Xn\Yo = yo)< Y K-kP\k\^A,yo) (B4) 

0<fe<n 

where we have used the notation 

R*k = Rk + R*k-i (B5) 

Summing both sides of eq.(B.3) over yo gives, with the 
same notations, 

1< E K-kS{k\TA) (B6) 

0<fe<n 

Applying the generating functions technique, we find 
then 

^ < R*{s)S{s\rA) (B7) 
where the generating function of i?^, eq.(B.5), is given 

by 

R*{s) = {l + s) Y ^2ns'" = (1 + s)-R(s) (B8) 

0<n<oo 

If now the A particle is fixed at the origin, eq.(B.l) 
becomes 

P{Yn = 0\Yo = yo)= Y P{yn=0\Yk = 0)P\k\0,yo) 

0<k<n 

= Y Rn-kP\k\0,yo) (B9) 

0<fe<n 

Summing both sides of this equation over yo and turning 
to the generating functions, we find, instead of eq.(B.7), 
the following equation 

^— = R{s)S{s\0) (BIO) 
1 — s 

Now, on comparing it with eq.(B.7), we infer that 

(l-s)S(s|r^)>^(l-s)5(s|0) (Bll) 

which implies that, asymptotically, if A moves, the reac- 
tion integral for a given trajectory Ta is not smaller than 
half of the reaction integral when A is immobile. 
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This unexpected conclusion requires some comments. 
First, it can be noticed that the equahty in eq.(B.ll) 
can be reahzed in a particular example. In fact, assume 
that B performs a classical Polya random walk on a d- 
dimensional lattice: at each step it can jump with an 
equal probability 1 / {2d) at one of the neighboring sites. 
Now, wc choose a special trajectory Ta for A, consisting 
of jumps from the origin to one of its nearest neighbors, 
1, and returns: thus, at each even time n, A is at 0, 
whereas A is at 1 at each odd time. 

For such a trajectory eq.(B.l) reads 

P{Yn=Xn\Yo=yo)= J2 K-kP\k\^A,yo) (B12) 



0<fe<n 



where 



r: 



n—k 



= Rn-k = P{Yn = 0|n = 0) = P(r„ = l\Yk = 1) 

if n — fc is even 

Rn-k+l = P{Yn = l|n = 0) = P(r„ = 0\Yk = 1) 

if n — fc is odd 



since 



R2n = P{Y2n = ^{Yq = 0) = {2d)-' X 
X ^P(F2n-l = y\Yo = 0) = P{Y2n-l = l\Yo = 0), (B13) 

where the prime designates that we sum only over the 
nearest to the origin sites. Consequently, we can write 



P°n-k — P-n-k + R n— fc+1 



(B14) 



Summing eq.(B.13) over yo we obtain for the generat- 
ing functions 



1 



= R°{s)S{s\Ta) 



(B15) 



with 



R°{8) = 



-R{s) 



(B16) 



Thus, for this special trajectory, we actually obtain 



(1 - .s)5(s|r^) = l{l-s)S{s\0) if s 



1 



(B17) 



which implies that the reaction is twice slower than for 
an immobile A. For instance, in three dimensions, the 
expression in eq.(B.16) is bounded when s ^ 1, which 
means that if A moves according to the previous rules, 
the reaction constant in case of a mobile A is half of the 
reaction constant for an immobile A. If an average is 
taken over trajectories, the Pascal principle can be valid 
or not, depending on the probability weight of the differ- 
ent trajectories. 

The physical reason of this conclusion may be under- 
stood if we consider the relative motion with respect to 



A. It is seen that at each integer time the relative dis- 
placement of B in the direction 0-1 can be or 2 at odd 
times, and or -2 at even times. On the other hand, the 
reaction integral at time n is related to the number of 
distinct sites visited by B up to time n. This number is 
clearly lower if A moves according to the foregoing rules, 
than if A is immobile, which explains that the reaction 
is slower in the first case. The same behavior can occur 
each time A and B are performing Polya walks with the 
same jump times. 

However, it should be pointed out that when A and B 
both perform Polya random walks with the same jump 
times, it may well happen that they exchange their po- 
sitions during simultaneous jumps. In this case, they do 
not react according to the rules we used here, but it can 
be relevant to adopt different rules, depending on the ac- 
tual phenomenon to be modelled. Then, the results could 
depend very much on these rules. 
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